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The zonal harmonic P« {6\ for large values of n and small values of 6, is 
approximately given by the formula 

P^(^) = Jo(a^), 
where a? = ?^ (^4- 1). 

This result is easily derived from a consideration of Legendre's equation, 

-p + cot ^ , -^+^(^ + 1)^ = 0, it = P^(5). (1) 

The first approximation of cot 6 is 1/5, and substituting this value in (1), 
the equation becomes 

of which a solution is Jo {ad), where a^ = n{n + l\ a result similar to the one 
discovered by Mehler,^ 

LtP„(cos5) = Jo(^). (3) 

4 -*■ CO \ ^/ 

The second solutions of (2), Bessel functions of the second kind or 
Neumann functions, do not appear in the expression for P«(5). 

Macdonaldf has given two formulae for Pm(5) in terms of Bessel functions. 



If ^ = (2#4"l)sin ~, 

P45) = JoH (4) 

* F. G. MeMer, * Math. Annalen,' vol. 5, pp. 140-144 (1872). 

t H. M. Macdonald, *Proc. London Math. Soc.,' vol. 31 (Old Series), p. 269 (1899) ; 
ihid,^ vol. 13 (New Series), pp. 220-1 (1914). 
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Bmaa 

P„(^) = Jo(i*)+sm2|||j3(a;)-J2(a;)+^Ji(a>)| 

12 0^ 8,':!:^^ J 

whilst Lord Eayleigh* obtained the result 

P» (^) = (1 + tI-J Jo («)+ 5^, Jo' (^). (6) 



f X^ 11^ 



12^^ 



6 a" 



in which = a^ and a^ = r^. (#4-1). 

Further terms in the coefficients of Jq(^) and Jo'(^) can be found without 
difficulty by successive approximations. Equation (6) may be written 



u = I\(^) 



Consequently 






and since cot 6 

equation (1) becomes 

(1% . 1 d% 



a {1 + ^3, '{ad), 

10 0^ 
3 



" 45 



= „(|+g)j.'(„.). 



Therefore 



#«^ ,1 d% 

-T — » -TT" *T* '?v 



l^^+iVo'(4 



or, the dependent variable u being changed to vJ^ (z), 






1 fz 



™ + -; 



(W 






cf02'5^\j0 '0/ a^\3 20^7 Jo' 

where the argument of the J functions is 0. 
Integrating (11), we find 



di) 






v_ 



dz = —-^ Ji^ 



dz zJq^ 

Now b^JoJo' 
It is easily seen that 



--a^f«*--I^h«')- 



2;' 



.2 



1 



0*JoJo »^ 



^ 



l(2Ji^-J32). 



f 






aW 






* Lord Eayleigh, ' Eoy. Soc. Proo,,' 1916, pp. 433-6. 
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Bmt 



and therefore 



izWdz = t^-l f .-/Ji^- J. )c?0, 



2 



f#Ji2^0=?^+'^ 



6 



3 



0*JlJ2<^^ 



Thus 



Hence 



I z'^JiJ2dz = I 0^ Ja . 2jJi ^2; = I z^J2 (^Js' + 2 Jg) dz 
= f (^^JsJb' + 2zW) dz^l U^J2^ = 






(13) 



^— ■'• f A ^^ J 2 . . ^ 



(2Ji2-J2^) k 



I? 



J. J "T" jCX 



^0 






6^^ J 



t^; 



120 #1 J? 



(14) 



From the recurrence formula, Jg^ may be expressed in terms of Jo and Ji : 

fT 2 
-—dz and 



^'?rf.+ 2^^'J^ 



J< 



Jo 



wZ, 



Integrating l ~~Y-^dz by parts, we get 
J Jo 



0^ji^ 7 r^^.^Ji^^ g/^Ji z^ 

^-rr4" «0 = 1 r:r^- dZ = SJ^f --i— ~ 

Jo a 



Jo=* 



Jo^ 



'j; 2/ 



Jo 2 . J Jo 4' 



and therefore 



a^J, 



*''^i'£^2+2|2^i&--^"''i 



28J1 



s* 



Jo^ 



Jo 



Jo 



Since B = 1 and A = 0, we find that 






z 



%a/ 



^ 



z^ 



Ji 






Z" 






3^4 



5a^ lOa^Jo \2 10 a^ 80 a' 



•^ 



and 



ti, = ( 1 4. 



z' 






z' 



12 a^ 160 a* 



+ ... 



z 



ma' 



:..JJof^) 



(15) 



(16) 



(17) 






... 1 Ji (2;) 



6a\ 10 oa^ ^ 



(18) 
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For the next approximation we obtain from (18) the value of dtc/dd as 
before. In this case 



When this is substituted, we find eventually the value of dv/dz, which now 
contains the term Iz^JqM^ as well as integrals having the product JoJo' in the 
integrand. But the former can be expressed in terms of the latter. Thus, 
integrating by parts. 



or, generally, 



1 






(19) 



Consequently, the problem of finding dv/dz depends upon the determination 
of the integral ^z^^JoJidz. 
Now 

j ^2»^ J^_ ^J^dz = I ^2»»- 1 J^ . zJ^^ idz 



\ 



^2»-iJ» (sJ™' + mJ™) dz = fH^ »' 



2 



(20) 



Also 



But 






pdz = z' 



^2m 



Z" 



2 



Z T 9 



f #'^-lJ/6f0 = ^ Jp 

I , 2m ^ 



2m-2p^2pJ^_^J^dz 

J^^^(^ftl^pJ L2m-~lJ^2d^^ 
1 



?-' mtf p Cv(4 



^2m 

^ -r 2 



<v -r 

2m 



z^"^ ^ 2 P 



and therefore 



2m 



(m+p) Iz^^^^Jp^z 



- \z^'^Jp(^Jp'-Jp+ijdz 

^2m- r^^dz + i f z^'^'ipip^idz, 



.-,2m 






J» + 



P 



/O ^ pO in-{^\CvZ>j 



(21) 



iZjJj 



or, combining the results in (21) and (22), 






m+j? 



m4"|? 



f^ 



p I'+i 



( 
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Applying this formula, making^ = 1, 2, 3, ,.., in succession, we obtain 



f #«go Jig0 






2m 



2 (m + 1) 

^2- (m-l) 

"17" "1 T — ~—~ 

2 (m-+l) 



/y2m r 



z^'''JiJ2dz 






>> 






J 2 . ^ T 2 

tn +-77* *^' 2 



(m + 2) L2 '^ 



a^^JgJgrfs 



r^4:^ 



Finally, we reach the term |0^'''J„i-iJ»/fe. Thus 



f 



P'^JaJidz 



2m 



J^2 



2(m-~l) j^2 ; 3(m->l)(m-2) j„2_ ' 
.m+1 (m4-l)(m+2) (m + l)(m4-2)(m+3). 



02»^ "" ^ = «*"-! (r+ l)(m- ~l)(m—2) ... (m—r) -y^ 



m + 1 L 



Ji^+ScosrTT' 



r = l 



(m+2)(w + 3) ... (m + r + 1) 



r+l 



. ( «uO^ 



Similarly, 



A<23» r 



m+p 



J? J/ + S COS rTT 



r = l 



(r + j9) (m —J?) (m —p — 1) . . . (m — p -^ r + 1) ^^ 



(m+J? + l)(^+J> + 2) ... {m-hp+r) 



r+p 



(26) 



f' 



3m- 1 J 2 



Z 



3m 



m +_p 



p 



■J 



1 r = m--|) 

2 m-p ,,ti 



• ' ~~ TT—. ^ rr:: —z ' ' — :: ' ti f^p 



(27) 



(m +p + 1) (m 4-p + 2) . . . (^t+j> + '?^) 
In the final step of each approximation, to obtain the value of t? it is 

2m+lT 2 



necessary to evaluate integrals of the type 
integrals occur under the form 



Jo tvZ an 1,1 I '^ 
0^ ' j Jo 



- clz and -^ — dz. These 



i 






&4-2m 



J 



dz. 



IIV ■ " • " v/j JL, ^,. . * « 







Integrating the first of these integrals by parts, and knowing that 



f 



"J? 



dz 



z3i z^ 



we obtain the result 

'^2»i+lJ 2 ^2ot+1J ^2m+2 



f 



0^0 



2m 



#»* Ji , 



(28) 



Jo "0 Jwi-f-J JJ0 

After further approximations, the expression for P»(^) in terms of J0(a^) 
and Ji {ad) is as follows : — 



P„W = («o-§ + 5-|+...)jo(«^) 



0_ 

6a 



(«i-§+g-|+-)ji(«^). 



(29) 
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where 



^0 



1 + 



2 01 



61 ^« 1261 5« 






79^10 



12 160 120960 29030400 20275000 



' "T" e • . 5 



'^*' 60 420 60480 
"^^ 315 560 ' 



S(i = 



2 






, , 6'2 31^* , 176)8 ^ 



9689^8 



10 3360 18900 106444800 



+ 



» 3 



Q _1, 461^31^*, 
^^ 5 105 6040 

~, — ^ _i. ^^ -i_ 

^^ 105 35 ' 



Si 



35 



To calculate P^ {6) to six places of decimals, when n is equal to or greater 
than 20, only a^, ^o, ax, and ^i need to be tabulated. The following tables 
have been constructed, giving these coefficients for each angle of the 
quadrant :— - 



r. 


«0' 


^0- 


%. 


^1. 





1 -0000000 


-00000 


1 -00000 


0-2000 


1 


1 -0000264 


0-00000 


1-00003 


0-2000 


2 


1 -0001016 


-00002 


1 -00012 


0-2001 


3 


1 -0002285 


-00004 


1 -00027 


0-2001 


4 


1 -m^mm 


0*00008 


1 -00049 


-2002 


5 


1 -0006350 


-00012 


1 -00076 


-2003 


6 


1 -0009146 


0*00018 


1-00110 


-2004 


7 


1 -0012453 


-00024 


1 -00149 


0-2006 


8 


1 '0016270 


-00032 


1 -00196 


-2008 


9 


1 -0020600 


0-00041 


1 -00247 


-2010 


10 


1 -0025443 


-00051 


1 -00305 


-2012 


11 . 


1 -0030801 


-00062 


1 -00369 


-2014 


12 


1-0036675 


-00074 


1 -00440 


0-2017 


13 


1 -0043067 


0*00086 


1 -00517 


-2020 


14 


1 -0049978 


0-00100 


1 -00601 


-2023 


15 


1 -0057411 


-00116 


1 -00690 


-2026 


16 


1 -0065368 


0-00131 


1 -00786 


-2030 


17 


1 -0073850 


0-00148 


1 -00888 


-2034 


18 


1 -0082860 


0-00166 


1 -00996 


0-2038 


19 


1 -0092401 


'00186 


1 -01111 


-2042 


20 


1 -0102476 


-00206 


1 -01232 


-2047 


21 


1 -0113087 


-00228 


1 -01360 


-2052 


22 


1 -0124237 


-00250 


1 -01496 


-2057 


23 


1 -0135930 


'00275 


1 -01636 


-2063 


24 


1 -0148168 


-00300 


1-01784 


-2069 


25 


1 -0160956 


-00326 


1 -01938 


0-2076 


26 


1 -0174297 


-00353 


1 -02099 


-2081 


27 


1 -0188194 


-00882 


1 -02267 


-2088 


28 


1-0202662 


0*00412 


1 -02442 


-2095 


29 


1 -0217675 


-00443 


1-02624 


-2102 


30 


1 -0233267 


-00475 


1 -02813 


-2109 


31 


1 -0249433 


0-00509 


1 -03009 


0-2117 


32 


1 -0266179 


-00544 


1 03212 


-2125 


33 


1 -0283508 


-00580 


1 08423 


-2133 
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e\ 


%• 


^0. 


ai. 


^1. 


34 


1 -0301425 


-00618 


1 -03640 


0-2142 


35 


1 -0319937 


-00657 


1 -03865 


-2151 


36 


1 -0339049 


-00697 


1 -04097 


-2160 


37 


1 -0358766 


-00739 


1 -04337 


-2170 


38 


1 -0379095 


-00782 


1 -04585 


-2180 


39 


1 -0400042 


0-00827 


1 -04841 


-2191 


40 


1 '0421613 


-00873 


1 -06104 


0-2201 


41 


1 -0443815 


-00920 


1 -06375 


-2212 


42 


1 -0466655 


-00969 


1 •05655 


-2224 


43 


1 -0490140 


-01020 


1 ^05942 


-2236 


44 


1 -0514277 


-01073 


1 -06238 


-2248 


46 


1 -0539074 


-01127 


1 -06542 


0*2261 


46 


1 -0564539 


-01183 


1 -06855 


-2274 


47 


1 -0590681 


-01240 


1 -07176 


-2287 


48 


1 -0617507 


0-01299 


1 -07606 


0-2301 


49 


1 -0645027 


0-01360 


1 -07846 


-2316 


50 


1-0673250 


0-01423 


1 -08193 


-2330 


51 


1 -0702185 


-01488 


1 -08551 


0-2346 


52 


1 -0731842 


-01654 


1 -08918 


-2361 


53 


1 -0762231 


-01623 


1 -09294 


-2377 


54 


1 -0793362 


0-01693 


1 -09680 


-2394 


65 


1 -0825246 


-01766 


1 -10076 


0-2411 


56 


1 -0857894 


-01841 


1 -10482 


-2429 


57 


1 -0891317 


0-01918 


1 -10898 


0-2448 


58 


1 -0925528 


-01997 


1 -11324 


0-2466 


59 


1 -0960538 


0-02078 


1-11761 


-2486 


60 


1 -0996361 


-02162 


1-12209 


0-2506 


61 


1 -1033009 


-02248 


1 -12668 


-2526 


62 


1 -1070496 


-02337 


1 -13138 


-2548 


63 


1 '1108835 


-02428 


1 -13620 


-2569 


64 


1 -1148042 


-02622 


1 -14114 


0-2592 


65 


1 -1188132 


-02619 


1 -14619 


-^616 


m 


1 -1229117 


0*02719 


1 -15136 


-2639 


67 


1 -1271016 


-02821 


1 -15667 


0-2664 


68 


1 -1313846 


-02927 


1-16209 


0-2690 


69 


1-1857621 


0-03036 


1 -16765 


-2716 


70 


1 -1402361 


-03147 


1 -17334 


-2744 


71 


1 -144808 


-0326 


1 -1792 


0-277 


72 


1 -149481 


^0338 


1 -1851 


-280 


73 


1 -154256 


-0360 


1 -1913 


-283 


74 


1 -159134 


-0363 


1 -1976 


0-286 


75 


1 -164119 


-0376 


1 -2039 


-289 


76 


1 -169212 


0-0389 


1 -2106 


0-293 


77 


1-174416 


-0403 


1 -2172 


0-296 


78 


1 -179733 


-0417 


1 -2241 


0-300 


79 


1 -185166 


0*0432 


1 -2312 


0-303 


80 


1-190716 


-0447 


1 -2384 


-307 


81 


1 -196386 


0-0462 


1 -2457 


-311 


82 


1 -202179 


0-0478 


1 -2533 


-315 


83 


1 -208098 


-0494 


1-2610 


-319 


84 


1 -214147 


-0611 


1-2690 


0-324 


85 


1 -220327 


-0629 


1-2771 


0-328 


86 


1 -226641 


-0647 


1 -2855 


0-333 


87 


1 -233094 


-0565 


1 -2940 


0-338 


88 


1 -239688 


-0584 


1-3028 


0-342 


89 


1 -246427 


-0604 


1 -3117 


0-348 


90 


1 -253314 


-0625 


1 -3209 


0-352 



The formulse and tables are of course especially useful when ^ is a small 
angle. They have been tested by computing P» {6) for comparatively small 
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values of n and large values of 9, Even in these extreme cases the results 
are correct to six places of decimals. 

For the purposes of comparison, a short table of V^{6) has been computed, 
giving in the second column the values of this function when the angle 6 is 
large, both 70 and 71 being required in these calculations. In the third 
column are the values calculated by Tallquist.* 





Ps {6). 




e. 


Calculated from (29). 


Tallquist's values. 


60= 


-0-073638 (9) 


-0-0736389 


65° 


0-241143 (9) 


0-2411438 


70° 


-0-278015 (2) 


-0-2780152 


75° 


-0-170219 (8) 


0-1702200 


80° 


+ 0-023307(8) 


+ 0-0233079 


85° 


+ 0-201746(2) 


+ 0-2017462 


90° 


+ 0-273437(6) 


+ 0-2734375 



In the construction of this table Jo (ad) and Ji(a^) were calculated from 
MeisseFs tables. When aO exceeds 15'5, Jo and Ji can be found to seven 
places of decimals from formulsef derived from the addition theorem of these 
functions. 

* Tallquist, * Acta Soc. Sc. Feniiicse' (Helsingfors), vol. 33, No. 4, pp. 1-8. Tables of 
P^ (cos B) for 71 ~ 1, 2, 3 ... 8 to ten places of decimals for each degree of the quadrant ; 
A. Lodge, Appendix to Lord Eayleigh's paper, " On the Acoustic Shadow of a Sphere/' 
Tables of Zonal Harmonics for 5° intervals in the value of ^ from 0° to 90^ and 
Qi = 1, 2, 3 ... 20 to seven places of decimals, 'Phil. Trans.,' A, vol. 203, pp. 87-110 
(1904). 

t "The Addition Theorem of Bessel Functions of Zero and Unit Orders, '^ ^Phil. Mag.,^ 
vol. 36, September, 1918. 



